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Scaling concepts and renormalization group (RG) methods are applied to a simple linear model of human posture control 
consisting of a trembling or quivering string subject to damping and restoring forces. The string is driven by uncorrelated 
white Gaussian noise intended to model the corrections of the physiological control system. We find that adding a weak 
quadratic nonlinearity to the posture control model opens up a rich and complicated phase space (representing the dynamics) 
with various non-trivial fixed points and basins of attraction. The transition from diffusive to saturated regimes of the linear 
model is understood as a crossover phenomenon, and the robustness of the linear model with respect to weak non-linearities 
is confirmed. Correlations in posture fluctuations are obtained in both the time and space domain. There is an attractive 
fixed point identified with falling. The scaling of the correlations in the front-back displacement, which can be measured in the 
laboratory, is predicted for both the large-separation (along the string) and long-time regimes of posture control. 
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I. INTRODUCTION 



A wide variety of systems subject to noise, random forces and interactions can be studied in depth by means of non- 
equilibrium statistical mechanics. This holds true whether the system in question is fundamentally of a chemical, biological, or 
physical nature. When the phenomena under study admits mathematical modeling by means of stochastic partial difi^erential 
equations, many powerful techniques can be used to analyze the effects that noise, fluctuations and random disturbances have 
on the dynamics as one changes both the spatial and temporal resolution scales at which the system is observed. The possibility 
to be able to use such techniques becomes especially pressing given that many typical real systems of interest are characterized 
by having many degrees of freedom interacting nonlinearly, leading to the competition between different length and time scales, 
with all scales evolving in the presence of noise and subject to uncontrollable external effects and contingencies. One of these 
important techniques is provided by the renormalization group (RG), suitably extended to dynamical systems and systems 
out of equilibrium |l||2j. Some recent results of renormalization group analyses of the kind presented in this paper have been 
obtained for diverse phenomena ranging from stirred fluids Q and turbulence fl^ , to surface growth phenomena p|-px|] , flame 
front propagation |T^], and cosmological large-scale structure formation p^p^ 

Fluctuations and noise are known to be present in physiological systems as well. Recently, a simple continuum model of 
human posture control was proposed jfTj that captures the gross or coarse-grained features underlying the physical mechanisms 
and adjusts well to laboratory measurements of time- varying displacements of the front-to-back (anteroposterior) sway recorded 
for human subjects in an stationary upright stance Despite the fact that the actual human postural control system 

must undoubtedly be highly complex, the stochastic model introduced by Chow and Collins ||f^ is described by a linear, and 
hence, exactly solvable, stochastic differential equation in one spatial dimension (in the following, temporal derivatives are 
denoted by an overdot, and the primes stand for spatial derivatives): 

j3y + y - vy" + ay ^ r){x,t), (f) 

where y{x,t) denotes the time-dependent front-back displacement measured with respect to the vertical located at j/ = (here 
we take the x-axis to denote the vertical axis as z will be reserved for the dynamical exponent which we introduce and calculate 
below). This is of course immediately recognized as a one-dimensional wave equation subject to friction (y), a linear restoring 
or pinning force (ay) and a stochastic or random force (rj). The onset-of-damping time scale is set by (3, the pinning time 
scale by and u is an effective string tension parameter. (This continuum equation describes the motion of a quivering 
damped string with a linear restoring force. A discrete version of this model would involve a chain of coupled random walkers, 
or monomers, and in this guise can be denoted as a pinned polymer [|I^). The stochastic noise is taken to be Gaussian with 
zero mean. For subsequent calculational purposes, it is convenient to characterize the noise spectrum directly in Fourier space 

{k,ujy. 

{ri{k,L0)'n{k',Lu')) = r{k,Lj){2TTf5{k + k')3{L0 + uj'), (2) 

where the angular brackets denote averaging with respect to the noise, and the noise spectrum function r(fe, uj) may contain 
in general both short range and long range correlations in space and /or time. We consider white noise here, so the spectral 
function is proportional to a constant noise amplitude: F — 2A. In IfTl, noise exhibiting short-time temporal correlations but 
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uncorrelated (or white) for long time scales was used in computing tiie correlation and response functions directly from the 
Fourier transform of the wave equation Eq.(|l]) where the former were fit to laboratory posture data allowing a phenomenological 
determination of /J,a~^, an effective noise amplitude and a parameter characterizing the short-time noise correlations. 

In 10 it is claimed that nonlinearities are not necessary to explain the posture data of healthy standing individuals, though 
it was also recognized that this may not be the case for subjects with balance disorders (nor even for normal individuals subject 
to a sufficiently large perturbing "kick"). These effects clearly lie outside the scope of the linear model Eq.(^) as do posture 
displacements outside the sway envelope or cone of stability [ poy2l| , and there is clear motivation for extending that model to 
include weak nonlinear terms. The minimal nonlinear term one can include in Eq.(|l|) is of the form ^ y^, which also serves to 
break the y — > —y symmetry in Eq.(^. Including this quadratic term makes good sense from the physiological point-of-view 
since real anteroposterior motion is intrinsically asymmetric [Q. Physically, this amounts to having a "pinning force" that 
varies with the amplitude of the horizontal displacement y. We will see, moreover, that such a term is needed in order to 
account for falling. In this paper we analyze the importance and impact of weak nonlinearities and to check the robustness 
of the linear model by means of a dynamical renormalization group analysis. We therefore consider the following nonlinear 
stochastic wave equation given by 

Py + y-uy" + ay + ey'^ ^r]{x,t), (3) 

where e > is the strength of the quadratic nonlinearity and we take the noise spectrum to be white (i.e., uncorrelated) at 
large scales: that is r(fc, ui) = 2A, where the noise amplitude is denoted by A, and serves also as a loop-counting parameter 
[pi] (an expansion of the solutions of Eq.(^ in loop diagrams is a convenient and powerful way to organize the calculation.). 

We are interested in the correlations in the solution(s) of Eq.(^ in the so-called hydrodynamic limit corresponding to large 
spatial separations (along the vertical axis) and long time intervals: \x — x'\ —> oo and \t — t'\ — > oo. This will tell us how 
the posture fluctuations are correlated along the length of the body at any given instant and how they are correlated in time 
at any given point on the body. In terms of Fourier variables (momentum k and frequency u>) this limit is taken by letting 
(fc, w) — > 0. The scaling information and universality class of this nonlinear wave equation is contained in two critical exponents: 
the dynamic exponent z and the "roughness" exponent x- These exponents are first obtained via a simple scale transformation 
of the stochastic equation of motion, and forms part of the full RG transformation (a course- graining or thinning-out of the 
degrees of freedom, followed by a re-scaling) A change of space and time scales x — > sx, t — > s^t is accompanied by a 
corresponding change of scale in the displacement field variable y s^y. Under this scaling the stochastic equation of motion 
Eq.(^ transforms according to 

— 2/3.. I . 2 — 2 II , z I z-^y 2 ifz — 1) — Y / i\ 

s py + y ~ s vy + s ay + s ^ey = s 2 ^ ' ^rj, (4) 

where we have used the the noise two-point correlation function Eq.(^ to determine the scaling of the noise source. Under this 
transformation, the individual parameters appearing in Eq.(^ therefore scale as 

/3^s-^/3, 

z-2 
V ^ S V, 

a s^a, 

^^s-2x-i_4, (5) 

At a fixed point of this scaling, the model parameters and the field y no longer change under a re-scaling for certain specific 
values of the exponents z and x- The model parameters approach their fixed-point values, P ^ P* ,v ^ v* , etc., and this fact 
gives rise to a corresponding fixed-point equation of motion, which is Eq.(^ written in terms of the fixed point parameters. 
Thus, each fixed point corresponds to a distinct dynamics governing the long time and large distance (s 00) behavior of the 
model. The dynamical phase space is thus divided or "partioned" into various domains or basins of attraction (or repulsion), 
each domain associated with a given fixed point. We can use this fixed point information to predict the asymptotic scaling of 
the displacement correlation function in both the temporal and spatial domains. This will be one of the main objectives of this 
paper. 

Independently of the RG, and in preparation for the results to be obtained, it is useful to derive the general scaling form 
of the correlation function of transverse displacements. Under a global scale transformation, the displacement field transforms 
according to 

y{sx,s''t) = s^y'{x,t), (6) 

which merely states that under a space and time rescaling, the displacement field can, and generally does, transform into a 
distinct (hence the prime) functional form y' , apart from picking up an overall factor. Thus, a change of scale will generally 
change the function itself, unless one is in the scaling or power law regime. When the system is known to be in a scaling regime, 
then in fact y = y' and from Eq.(^), the auto-correlation function therefore scales as 
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{y{x, t)y{0, 0)) = s-^"" {y{ax, sH)y{Q, 0)), 

= .^x^(±), (7) 

where (without loss of generahty) we have chosen s ~ and 5' is a (dynamic) scahng function, which itself exhibits power 
law behavior for asymptotic limits of its argument: 

^ / X r A for u ^ 0, 

* « 2x ' (8) 

^ ' Bu ~ tor u —* oo, 

for constants A, B. The dynamic exponent z describes the scaling of relaxation times with length and x is the "roughness" 
exponent of the string or polymer. Thus, knowledge of these two exponents is all that is required to determine the explicit 
scaling of the correlation function within each dynamic phase of the model. For the linear model (i.e., e = 0) these exponents 
can be exactly determined with little effort, and there are just two solutions (for white noise). For one, the wave equation for 
non- vanishing pinning force (q 7^ 0) is made scale-invariant with the choice z = and x = This exponent solution can 
be read off directly from Eq.(p|) taking into account the fact that the noise amplitude is constant and hence, non- vanishing on 
all scales (.4 > 0). This immediately yields the exponent identity z = 2x + 1. For a finite fixed value of the pinning force, the 
only possibility is to take z = 0, since a positive z > yields an asymptotically divergent a, while 2; < would yield instead an 
asymptotically vanishing value. We see that /3* = /3 is finite and the diffusion constant vanishes, u* = 0. In this phase then, 
there is no diffusion. This exponent pair corresponds to the experimentally observed scaling regime denoted as "bounded" or 
"saturated", and holds for the very latest times when the pinning force has had time to correct for posture excursions from 
the vertical and aligns the body in an upright stance In the earlier "diffusive" scaling regime, the pinning force has not 
had sufficient time to act and is negligible, i.e., a ~ 0, and there is another exact exponent solution given by z = 2 and X ~ ^) 
indicating that in this parameter regime, the model belongs to the same universality class as the one-dimensional Edwards- 
Wilkinson (EW) model [^. At this fixed point, we have a finite diffusion u* = ly and /9* — 0. In this phase, there is no wave 
propagation, since the second derivative in time is absent. Note that these simple scaling solutions have been obtained from 
applying naive scaling arguments to the linear equation. However, as soon as the nonlinearity e is turned on, and no matter how 
weak, other nontrivial exponent solutions arise for which the nonlinearity can become relevant. The naive scaling arguments 
are insufficient for obtaining the scaling exponents in the fully nonlinear model. The RG allows one to calculate z and x in the 
combined presence of fiuctuations and nonlinearities and to calculate the exact asymptotic scaling of the correlation function 
Eq.(0) in all the basins of attraction. 

The rest of this paper is organized as follows. In Section II we make use of a dynamic functional formalism for the perturbative 
calculation of solutions to Eq.(^ based on the Martin-Siggia-Rose (MSR) Lagrangian. The bare correlation function, response 
function, noise spectral function and bare interaction vertex function are identified, their corresponding Feynman diagrams are 
introduced and calculated and provide the basic elements of a systematic and controlled loop expansion for the one-particle- 
irreducible (IPI) diagrams which we then use for extracting the one-loop RG equations in the low-energy regime. 

In Section III we exhibit the set of nonlinear differential RG equations for the dimensionfuU parameters appearing in Eq.(|^. 
For white noise these involve five equations: one associated with each independent parameter appearing in the equation of 
motion. We then identify a convenient set of three dimensionless couplings in terms of which these RG equations can be 
expressed. The RG flow is therefore represented in a three-dimensional dimensionless parameter space and we solve for all the 
one- loop fixed points in terms of this reduced set. In this way we find a total of two trivial fixed lines and four non-trivial 
fixed points. Linearization of the RG about each fixed point (or line) reveals the nature of the fixed point, in the dynamical 
systems sense (whether the fixed point is a source, a sink, a limit cycle, a spiral, a saddle point, etc.) and yields linear stability 
information which we quote in terms of the eigenvalues and eigenvectors of the linearized RG. The two lines of fixed points 
correspond to the diffusive and bounded phases of the strictly linear model, and are present in the non-linear model for all 
values of a and /3. Of the four non-trivial fixed points, one is a stable spiral which represents the "falling" phase. The other 
three are saddle points which seem to have rather little infiuence on the long-range and long-time dynamics, however. 

Substitution of the fixed points back into the original set of RGs yields the values of the critical exponents z, x for each 
fixed point and hence determines the exact asymptotic scaling properties of the correlations Eq.(0) in the basin of attraction 
(or repulsion) of each fixed point, which is presented in Section IV. The numerical analysis of the fixed points is then repeated 
using an alternative set of three dimensionless parameters suitable for investigating the small a limit (corresponding to the 
diffusive regime). The use of this second set of parameters in conjunction with the first is necessary in order to completely 
cover the entire model parameter space. 

The detailed structure of the nonlinear RG fiow is revealed by plotting the fixed points and mapping the numerically computed 
(and normalized) vector field of the non-linear RG fiow in the neighborhoods of all the points in Section V. Many important 
aspects of the morphology of the dynamic phase space are qualitatively revealed and allow conclusions to be drawn regarding 
the impact of the weak nonlinearity. This provides revealing information regarding the shape and structure of the basins of 
attraction and complements the analytic analysis. Summary and conclusions are drawn in Section VI. 

A number of explicit analytic calculational details needed for the derivation of the RG equations are relegated to the 
Appendices. The complete calculations leading to the one-loop response function are presented in Appendix A, and similar 
calculations for the noise spectral function and vertex renormalizations are given in Appendix B and C, respectively. 
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II. DYNAMIC FUNCTIONAL FORMALISM 



In this Section we make use of a functional integral representation of non-equilibrium stochastic dynamics. This leads to the 
efficient identification and extraction of the calculational elements (and Feynman rules) needed for the perturbative calculation 
of the solutions of any stochastic partial differential equation. It is well known how to map stochastic ordinary or partial 
differential equations with additive noise into equivalent generating functionals | p^j23| . Essentially, there are two formally 
distinct but physically equivalent routes one may follow, an option one has at least in the case of Gaussian noise. In the 
Martin-Siggia-Rose (MSR) formalism p^-p^ one introduces a fictitious conjugate field (call it y) with its own source term. 
The equation of motion, in this case given by Eq.(^, is imposed as a constraint on the dynamic functional, and is realized 
linearly. In the minimal formalism, no conjugate field is introduced, leading to a nonlinear realization of the constraint | p3| , p2[ . 
For Gaussian noise in the minimal formalism, the constraint appears quadratically in the argument of the functional, while for 
non-Gaussian noise, it is inherits whatever nonlinearities are present in the noise probability distribution function itself [ p2[ . 
Here we develop the calculation following the MSR approach, since this leads to a simpler structure for the associated Feynman 
diagrams, and our immediate aim is to obtain a perturbation expansion which can be set up, organized and calculated in terms 
of a few elementary diagrams or graphs. 

The MSR dynamic generating functional corresponding to Eq.(^) is given by (taking a translationally invariant noise spec- 
trum) 

Z[J, ,]]= j [dy] [dy] exp - i J dxdt yVy + i j dxdt y {Py + y - uy" + ay + ey"^] + j dxdt [yj + yj]^, (9) 

where y denotes the conjugate field and J, J are arbitrary sources for y and y, respectively. The noise r] has been integrated 
out exactly, and appears in this functional only through its two-point or correlation function F. There is also in principle 
a certain Jacobian determinant factor in passing from Eq.(^ to Eq.(|^, but it can be shown on general grounds to be a 
constant, and hence irrelevant for computing normalized correlation functions (see, e.g., ^,^). The noise spectrum F as 
written here is understood to be given in terms of x,t. All the dynamic and fluctuation information contained in Eq.(^ is 
also contained in Z, which is an alternative representation of the dynamics. In preparation for the RG transformation, which 
is most straightforwardly implemented in the Fourier domain, we cast this functional in terms of momentum and frequency 
variables from the outset. 

To this end, we introduce Fourier transforms for the physical and conjugate fields and the noise spectrum, i.e., 

r§ r^y(fc,..)e^("\ (10) 

<J J — oo 

where, in a mild abuse of notation, we distinguish the functions from their Fourier transforms only through their arguments; 
this however avoids a clutter of notation later on. Note we implicitly cut-off the momentum integration at the scale A = ^ 
where a plays the role of a minimum distance of spatial resolution or lattice spacing. The cut-off sy mbo l on the integral (>) 
means one is to integrate over all momenta in a "shell" such that A/s < |A;| < A where s > 1 (see Eq. |A6|) . The cut-off defines 
the spatial scale above which it makes sense to use continuum equations for modeling. The quadratic or Gaussian part of the 
functional (i.e, e = 0) can be exactly computed and serves as the starting point for aperturbative expansion which we will use 
to calculate the RGE's associated with the nonlinear stochastic wave equation Eq.(y). From standard Gaussian integrations 
[ pilJi^ , we have that (up to an overall irrelevant constant prefactor) 

^rr ri f ^ dk T dio f 1 J(fc, ^)F(-fc, -c^) J(-fc, -^) , . J {k,uj)J{-k, -uj) \\ 

Zo[J,J|-exp(^y — y^^— I ^2 + [^fe2 _ + «p +'iu;-l3Lo2 + uk^ + aj)- ^^^^ 

From Eq.(^), it is clear that all noise- averages of arbitrary products of physical fields and conjugate fields at distinct points 
and times are obtained from the appropriate functional derivatives of ln(Z) with respect to the source terms J, J, taking the 
sources to zero at the end of the calculation. We may thus obtain the "bare" or zeroth order auto-correlation and response 
functions directly in Fourier space (the zero (0) subscript denotes the zero-coupling limit e = 0) as follows: 



{y{pi,uji)y{p2,uJ2))o 



(27r)^ S^Zo[J,J] 

Zo[J,J] SJ{pi,UJl)5J{p2,UJ2y^=^=° 
= (2tv) d[pi +p2)d(LJl +LJ2) 2 2 12' 

= {2TvfS{pi + P2)S(lOi + LU2)C(pi,iUiy, (12) 



for the auto-correlation function and 
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/ ( ^v ^\ (27r)" S^Zo[J,J] 

Zo[J, JJ ()J(pi,tJl)dJ(p2,l^2) 
_ .{2TvfS{pi +P2)S{UJ1 +LJ2) 



iuji — I3uj\ + vp\ + a ' 

= {2Tvf5{pi + P2)S{lOi + CJ2)iA{pi,UJl), (13) 

for the response function. The final equations in each case define the corresponding reduced auto-correlation and response 
functions: C and A, respectively. The delta function factors reflect the fact that momentum and energy (or frequency) 
are conserved as they "flow" through the correlation and response functions. It will be noted that the reduced response 
and correlation functions are related via C = A'FA — r|Ap, and this relation will be used to obtain the noise amplitude 
renormalization from knowledge of the correlation function. They also automatically satisfy the fluctuation-dissipation theorem: 
C{p, uj) = ^ ImA(p, —Lu). The remaining element from which the complete perturbative expansion of the dynamic functional is 
developed is provided by the vertex function, which necessarily involves the trilinear interaction term. The complete dynamic 
functional Eq.(^ can be written as 

Z[J, J] = exp (^[A, Zo[J, J]^J2 ^ ^"t-^' -^l' (14) 



n=0 



where the exponential of the following interaction operator acts on the Gaussian part (Zq) of the generating functional: 

Si[4-r, ^] = / I TT \ (2vr)'<5(fci + k2 + k2)5{iOi +i02+ i03) {2nf—. , (15) 

'SJ'Sj' J (2^)' J ' SJ{kl,UJl)SJik2,u;2)SJ{k3,iU3) 

and merely reflects the fact that the non-linear interaction term for the MSR functional is a cufeic-interaction involving one 
conjugate field and two factors of the physical field: ~ ieyy^- The delta functions are a consequence of the energy and 
momentum conservation at the vertex, and the strength of the vertex is given by the parameter e, which in fact defines the 
reduced vertex. We complete the specification of the elements encountered in a perturbative expansion of Z by specifying the 
vertex function, which we read off simply by inspection of Eq.([l5|): 

e(fci, fc2, fcs; oJi, t^s) = (27r)^5(fci + k2 + k2)S(uJi + UJ2 + i-L!3)ie. (16) 

Each of these calculational elements, C, A, F, e, can be represented by a simple diagram as depicted in Figure (Q), and then 
used to systematically construct a perturbation expansion of the dynamic functional, and its Legendre transform (or effective 
action) which is the generator of the one-particle irreducible (IPI) diagrams. (Note: we remove the factor of i appearing in 
the bare vertex (^5|) and reponse functions ( |l^ ) by re-defining the conjugate field y ^ iy m the functional integral.) The 
Feynman diagram for the bare response function A is shown in Figure lA, which reflects the fact that the response function is 
constructed from the mixed product of the conjugate field (wiggly line segment) times the physical field (straight line segment). 
The arrow indicates the direction that momentum k and frequency uo fiow through this diagram (from the conjugate to the 
physical field) and is the convention we adopt here. It is important to keep this in mind since A(fc, w) 7^ A(— fc, —uj), as can be 
seen by inspection of (^^. The bare correlation function is shown in Figure IB. As this is an even function in both momentum 
and frequency, we need not indicate a flow direction for this function. Finally, the bare trilinear vertex is shown in Figure IC. 
This involves one conjugate fleld and two physical fields. The bare interaction is represented by a small open circle. The arrows 
indicate the fiow direction for the two physical fields and is just a convention; the crucial point is that momentum and frequency 
are conserved at each vertex, as per (^|). These are used to build up the IPI diagrams, which are those diagrams that cannot 
be broken down into disconnected sub-diagrams by cutting an internal line (see examples below). The class of IPI diagrams are 
of fundamental importance in any perturbative scheme based on diagrams since all other diagrams can be unambiguously built 
up from these primitive "building blocks" . In Figures 

%i (§),(§), we depict the IPI one-loop diagrams contributing to the noise 
spectrum, response function and vertex function. These diagrams are built up from the elementary vertex and response and 
correlation functions. Note that the bare correlation function is itself a composite function built up from the noise spectrum 
and response functions and is not IPI. This fact is also reflected mathematically in the fact that the bare correlation, noise and 
response functions obey a fluctuation-dissipation theorem. The complete mathematical transcription of the diagrams is carried 
out in the Appendices. 

Finally, we note that each loop diagram is multiplied by an overall symmetry factor. This factor receives a contribution from 
the factorial coming from expanding out the exponential prefactor in ( p^ ) to a given order in the coupling. For example, a 
graph with n bare interaction vertices yields a factor of ^ . The other contribution comes from simply counting the number of 
distinct ways in which the given diagram can be built out of elementary Feynman graphs (Fig. lA-C) keeping the topological 
structure fixed Multiplying these two numbers together yields the final net symmetry factor. For the response, noise and 
vertex one-loop diagrams, these turn out to be 5*^ = 4, Sn ~ 2 and Sv ~ 4, respectively. 
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FIG. 1. A: The bare response function, B: The bare correlation function, C: The bare vertex function. 



/VM>^ = Sr 

FIG. 2. The one-loop corrected response function. 




FIG. 3. The one-loop corrected vertex function. 

o = o + Sn { AA/K^:>vv^ } 

FIG. 4. The one-loop corrected noise spectral function. 





III. DYNAMIC RENORMALIZATION GROUP EQUATIONS AND THEIR FIXED POINTS 



By means of the diagrammatic expansion, we calculate the one-loop corrections to the response function A, the noise spectral 
function F and the nonlinear vertex e and from them, obtain the associated one-loop corrections to the model parameters 
a, P, u, A and e in the large-distance and long-time limit. The complete one-loop expressions are derived and calculated explicitly 
in the Appendices A, B and C for the response function, noise spectral function and the interaction vertex, respectively. With 
these in hand, the next task is to carry out the renormalization group transformation, leading to the differential flow equations, 
whose steps we briefly review here [Q]. 

1. We first perform an infinitesimal Kadanoff or "block"-transformation: that is, we integrate over a thin momentum shell 
A/e* < |fc| < A where I— l~'r5,Q<5 << 1. This means we integrate over (small-scale) fluctuations characterized by 

having their momentum in this range: (/_^ + Physically, this step serves to thin out the degrees of freedom 

(coarse-graining) and reduces the spatial resolution of the system. Note that in one space dimension, the "momentum 
shell" reduces to two disjoint intervals. For d > 2 and higher, it is a thin spherical shell in momentum space. 

2. After performing this step, the resulting equations have a new, lower momentum cutoff of A/e*. This means we have 
changed the lattice constant of the system, so to restore it to its original size, we re-scale the momenta by putting 
k — > ke~^. This is identical to the scaling carried out earlier with s = . The parameters are rescaled (but here, 
differentially) as in , but with additional corrections coming from the momentum-shell integrations carried out in step 
(1). 



Applying these steps to the one-loop expressions (A8, B5, 04), we obtain the following set of differential renormalization group 
flow equations: 



da 
U 

d£ 
'dl 



^(-^-(rw^-^ + 



45 



2/(1 + ^)4)' 



3/ 



f{l + h)-l \ {l + h) 4(l + /i)2 4(l + /i)3 



de ( 



dv 
'dl 
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(l + /i)2 t(l + /l) 

2<7 



/ + 



( — 

\2(l 



3f 



4/(1 + h)-l 
f 



[2 



4/(1 + h)-l 
8hf , ( f 



[4/^(1 + /0 +3/ 



+ 



4/(1 + h)- 
1 



(i+ft) 

5/ 1 



(l + /i) 4(H-ft)2 4(l + /i)3 



{l + hY {l + hf 4(1 + ft) 



where the three dimensionless parameters g,h, f are defined by 



9 - 



e'AA 



h = 



f = a/3. 



(17) 



(18) 



This reduction from five dimensionfuU to three dimensionless parameters is a concrete realization of a more general result 
from dimensional analysis known as Buckingham's D-Theorem pol], which states that if one has m dimensionfuU variables in 
a theory involving n fundamental units (such as length, time, mass, etc.) then there exist m — n independent dimensionless 
groups or combinations of the m original quantities. In our case, m = 5, and n = 2, since all the parameters appearing in 
the stochastic equation (^ can be expressed in terms of two fundamental units: namely length and time. Thus, we expect to 
be able to write the five RG equations exclusively in terms of 5 — 2 = 3 independent dimensionless groups or combinations of 
the original dimensionfuU parameters. The fact that a cut-off is introduced from the RG transformation presents no problem 
since [A] — has units of inverse length. The number of fundamental units remains the same. Below we introduce a 
second, independent group of three dimensionless parameters which when taken with this first set, will suffice to cover the 
entire parameter space. 

The dimensionless coupling controlling the nonlinearity is g, and a non-zero value of this coupling is what drives the RG fiow 
away from the free, or Gaussian limit. Since g oc e^A, to have nontrivial fiow simultaneously requires both a non-vanishing 
nonlinearity e 7^ and fluctuations ^ > 0. Indeed, by setting g = in (^), we immediately recover the naive scaling laws 
given by (|), expressed here in differential form, with s = . Thus we see that the model parameters depend on scale in a 
complicated way when nonlinearities [g 7^ 0) are present. In terms of these dimensionless variables, the RG flow is given by 
the following set of three differential equations: 
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d£ 



-J1J = 9 



[15 + /(I + h) + 6ft]) , 



4f 

M 
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(l + ft)3 

g [ 
(l + ft)-*' 

-(l + 4/(l + ft)2 + 2f (l + ft)2). 



3/i + 2/i^ + /(-I +2h + 3h^ 



(19) 



2(l + ft)4 

These are obtained by differentiating ( |l8| ) and using (|l7|). Note the dependence on the two exponents has dropped out. We 
solve for the complete set of real fixed points whose coordinates are {g* ,h* , f*). These are labelled, collected and presented 
in Table I along with their associated critical exponents, infrared eigenvalues, eigenvectors and the nature of the fixed point. 
There is actually an entire line of attractive fixed points, denoted by LI, located at (0, 0, /) , where / is any real number. A 
fit of the linear model to the data yields /* — 0.05. As we will see, LI is associated with the final saturated phase for which 
v* — 0. This implies that h* — 0. In the linear model g = g* = 0, so the posture data picks out a distinguished point on the line 
LI, and we see the final phase of posture control is controlled by an infrared attractive fixed point at (0, 0, /*). The quadratic 
term gives rise to three saddle points, P3, P4, and P5, and one stable spiral, P2; their coordinates are listed in the second 
column of Table I. When we substitute these fixed points into the original set of RGB's (|l^), we obtain the corresponding pair 
of critical exponents: z = z{g* , h* , /*), x = x(S'*i /*) associated with each fixed point (or fixed line): 



2g* 



(1 - 3*[3 + /*(1 + h*) + 2fa*] + /i*[3 + 3fa* +h*' 



(20) 



(l + /i*)2 ^ 2(1- 

which we obtain from using the pair of equations for A and a in ( |l7| ) to solve for z,x- Any pair of equations depending on z 
and X taken from (|l^) can be used to solve for these exponents provided that neither equation in the pair trivially evaluates to 
zero at the fixed point whose exponents one wishes to compute (thus, for example, if /3* = , then we cannot use the equation 
dj3/dl = to solve for z, and so forth.). Taking this obvious restriction into account, we have checked that the exponents 
calculated taking all possible pairs of the RGB's agree. This serves as an important consistency check of the entire calculation. 
The exponents are listed in the third column of Table I. By linearizing the RG equations ( |l9| ) around each fixed point, we 
obtain additional information regarding the stability of the fixed point and can classify the infrared stability properties of the 
point. This information is helpful for visualizing the RG flow and interpreting the flow graphs in Section V. At each flxed 
point we substitute g = g* + Sg,h — h* + 5h, f = /* + 5f, into ( p^ ) retain all terms up to order 0{Sg,5h,Sf) and compute 
the eigenvectors and eigenvalues of this linear system. This information completely characterizes the RG flow in the linearized 
neighborhood of each flxed point. The eigenvalues and eigenvectors are listed in columns four and flve of Table I. The unit 
direction vectors in this coupling space have components g — (1, 0, 0), h — (0, 1, 0), / = (0, 0, 1). Note that the line LI has only 
two RG flow eigenvectors associated with it. From knowledge of the eigenvalues, we immediately predict the nature of the 
flxed point, and this is given in the last column of Table I. 



Fixed 
Point 


Position 


Critical exponents 


IR-eigenvalues 


IR-eig 


envectors 




Class 


LI 


(0,0,/) 


(0,-1/2) 


(-2,-1,0) 


vL = 
vL = 


(1,0,0) 
(0,1,0) 




Stable line 
attractor 


P2 


(-0.37,-2.15,-1.79) 


(-0.56,-0.41) 


(-6.79, -0.67 ±0.27i) 


"P2 = 
VP2 = 


(-0.17,-0.99,0.03) 
(0.57 + 0.111,0.72,- 
(0.57-0.1H,0.72, - 


-0.24 - 0.28i) 
-0.24 + 0.28i) 


Stable spiral 


P3 


(0.07,0,-0.29) 


(0.14,-0.41) 


(-1.78,1,-0.10) 


vpa = 

"PS = 
VP3 = 


(0.06,1.00,-0.04) 

(1.00,0.00,0.00) 

(0.00,0.00,1.00) 




Saddle point 


P4 


(0.08,0,-1.71) 


(0.15,-0.45) 


(-1.65, 1., 0.11) 


UP4 = 
fp4 = 
WP4 = 


(0.08,1.00,-0.04) 
(1.00,0.00,0.00) 
(-0.01,0.00, 1.00) 




Saddle point 


P5 


(2.86,-2.96,-0.07) 


(1.48,0.03) 


(-8.24,1.50,-1.28) 


VP5 = 
2 

VP5 = 

Vp5 — 


(0.21,0.98,0.01) 
(0.97,-0.23,0.01) 
(0.79,0.45, -0.41) 




Saddle point 



TABLE I. 
Fixed Points in terms of g, f,h. 
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Of the four non-trivial fixed points, three of them {P2, P3, P4) are well within the limits of perturbation theory as they all 
correspond to small values of the dimensionless coupling \g*\ < 1. By contrast, the coupling for the P5 fixed point turns out to 
be rather large: g* = 2.86 > 1, and it is likely this point is an artifact of perturbation theory. In any event, as can be seen from 
inspection of the graphs in Section V, P5 is "far" away from the remaining fixed points and has very little or no impact on the 
RG fiow in the small coupling region, which is the region safely explored by perturbation theory. Both P3 and P4 are near 
the point (0,0,/*) belonging to LI, but a study of the flow in the neighborhood of (0,0,/*) shows these to be unimportant 
and substantiates the claim made in that nonlinear effects are not needed to explain posture data of healthy individuals. 
Nevertheless, nonlinearities are needed to account for falling, and the stable spiral P2 is an attractor for "falling" , as we now 
show. For P2 the fixed coupling g* = —0.37 is negative. But g* < =^ a* < 0, since e > and .4 > 0, and a change in 
sign of a signals the transition from stability of the upright vertical stance j/ = to instability or falling : y > 0. This can be 
seen quite clearly by examining the properties of the effective mechanical potential U{y) associated with our equation. This 
potential is obtained by integrating the y-dependent force terms appearing in (p|) and yields 



(21) 



since F{y) = -—U'{y) is the deterministic (non-random) force acting on the string. A plot of U versus y is shown for a > and 
for Q < in Figure 5, where U{y) is plotted along the vertical and y along the horizontal axes, respectively. 



U(y) 



U(y) 



0.2 

0.1- 


-0.1 
-0.2 



0.1 



-0 . 1 
-0.2 
-0.3 
-0.4 









[a] ; 


2-10 1 2 

y 








[b] 



FIG. 5. Effective mechanical potential U{y) for [a] a > and [b] a < 0. U{y) has units of energy and y has units of length 



For Of > 0, the origin at i/ = is locally stable to small perturbations (to fall, one would have to be pushed backwards 
sufficiently hard so as to overcome the potential barrier at j/ = —1. However, when a < 0, the vertical stance becomes an 
unstable configuration, and a fall in either direction results (a clear mechanical example of symmetry breaking). This falling 
phenomena is beyond the scope of the linear model since a is always positive in the linear model. The flow graphs in Section 
V reveal that P2 lies in a domain or phase which is separated from the domain or phase of upright stance. Thus, healthy 
individuals are characterized by having their initial conditions of posture control the basin of attraction of LI. 

The above group of dimensionless couplings forms a useful set when v is small, and the limit as v tends to zero (corresponding 
to vanishing viscosity or diffusion) can be safely studied. The limit of small /3 can also be treated with this set. When a is 
small however, corresponding to the diffusive regime, a different set of couplings is required. To study the flow corresponding 
to a ^ (corresponding to a vanishing linear restoring force) we introduce the set of dimensionless variables 

G = 9/h'-^, H^h^' = -^, F = f = ap. (22) 
In terms of these new variables, the original set of RG flow equations (|l^) takes the following form: 
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da ( 



2G 



Ti 

dA 

dl 
df 
dl 
de 
dl 



G \ F 



+ ■ 



H 



{l + HY\H 2F{1 + H) 



A{z-2x-l + 



F 



+ 



(l + //)2 [H (1 + H) 



— = v[z-2 



G 



— =e 2 + X + 



(1 + ^^) 
4G 



(23) 



By setting G — 0, we once again recover the naive scaling laws given by (^). In terms of this group of three dimensionless 
couplings, the RGE flow is characterized by 



H{1 + HY 
G 



[4F(-2 - H + H^)+ 6H{2 + H) 



de + (i + i/)4 



dF 

He 



F 



[-2 -3H - 3H^ + F{-3 - 2H + H'' 
GH 2G FG 



2F{1 + HY H{1 + HY H{l + Hy 



(24) 



These result from differentiating (|22| ) and using (^3|). Once again, the dependence on the two exponents z,x drop out of the 
equations for the dimensionless variables. We solve for the complete set of real fixed points which we denote collectively by 
(G* , H* , F*) and are displayed in the second column of Table II. We note that we reproduce two of the fixed points obtained 
with the original variables, namely P2 and P5. Although their coordinates as expressed with the variables (G, H, F) are distinct 
from those in the (<?, h, /) system, their exponents, eigenvalues and eigenvectors are identical in both systems. Once again, we 
see that P2 is an attractor for falling: in these variables G* > z/* > so that H* < =^ a* < 0. 

We also discover a repulsive fixed line 11 located at (0,0, F) — (0,0,/). This fixed line was "missed" by the first group of 
parameters. The posture data in the diffusive phase selects a point on this line, namely (0, 0, 0.05). 

We substitute these fixed points into the original set of RGE's (^^ to obtain the corresponding pair of critical exponents: 
z = z{G* , H* , F*),x ~ x{G* ,H* , F*) at each fixed point: 



2G* 



H*{l + H' 



X = 



2G*(l + 3f/*) 



(25) 



These follow from taking the equations for e and a. Taking other pairs of RGE's yields expressions for the exponents that 
evaluate to the same numerical values at the fixed points. In analogy with the first set of RG equations, we calculate the 
eigenvalues and eigenvectors associated with the linearization of ( p^ ) about each fixed point (fixed line) and determine the 
nature of each point from the eigenvalues. This information is organized below in Table II. 



Fixed 
Point 


Position 

(G*,H*,F*) 


Critical exponents 

(2-X) 


IR-eigenvalues 


IR-eigenvectors 


Class 


11 


(0,0, F) 


(2,1/2) 


(5,2,0) 


Vi\ = (1,0,0) 
Vfi = (0,1,0) 


Unstable line 


P2 


(0.04,-0.46,-1.79) 


(-0.56,-0.41) 


(-6.79, -0.67 ±0.27?) 


V^2 = (0.16,-0.98,-0.13) 

via = (0.05 - 0.02j, 0.25 - 0.30i, 0.92) 

Via = (0.05 + 0.02j, 0.25 + 0.30i, 0.92) 


Stable spiral 


P5 


(-0.11,-0.34,-0.07) 


(1.48,0.03) 


(-8.24,1.50,-1.28) 


Vps = (-0.72, -0.69,0.09) 
Vis = (-0.39,0.88,0.26) 
vis = (0.19,0.12,0.97) 


Saddle point 



TABLE II. 
Fixed Points in terms of G, F, H . 
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IV. CORRELATION FUNCTION: SCALING PROPERTIES 



As discussed in the Introduction, the RG fixed point analysis can be used for predicting the asymptotic, large-distance and 
long-time limits of the two-point correlation function of transverse displacements from the vertical upright position: 

C{xi - X2,ti-t2) = {{y{xi,ti)-y{x2,t2)f), (26) 

which measures the fluctuations in the difference of transverse displacements at two different points along the body and/or at 
two different times. Because the model is translationally invariant, this function depends only on the differences xi — X2 and 
ti — t2 in space and in time. In the scaling regime, which holds when the system is in the vicinity of one of its fixed points, it 
is easy to derive the exact scaling behavior of which emerges in the large distance and long time limits. Put x = xi — X2 
and T — ti — t2. Consider correlations in the time-domain measured at the same point on the body, so that a; = 0. Then from 
y{x,T) = s^^y{sx, s't) we have that 

C(0, r) = s-^^'CiO, sV) ~ r^'^/" C(0, 1) = B t^^^\ (27) 

with B a constant, which follows from choosing s such that s^t = 1. Next, for correlations in the spatial domain, set t = 
and we obtain 

C{x,0) = s-''''Cisx,0) ~ x^^C{l,0) = Ax'''', (28) 

which follows from choosing s such that sx = 1; yl is a constant. These results are consistent with and imply the scaling limits 
of the scaling function appearing above in (|7|). 

Using the Tables |l[ 0, it is a simple matter to calculate the scaling of the correlation function about each fixed point. First, 
consider the limiting case represented by "turning" off the quadratic interaction term: e ^ 0. In this limit, we recover the 
linear model of Chow and Collins (^ and we reproduce the corresponding scaling behavior(s) of the correlation function. In 
this limit (and for uncorrelated noise) there are then only the two trivial fixed points (actually, fixed lines) II, LI (see Table Wi , 
corresponding to the small a and small v limits, respectively. Reading off the critical exponents from Table ^ and using ( |2^ ) 
we confirm that in the neighborhood of 11 the correlation function scales as C(r) = Bt^^^. (In the temporal correlations 
are parametrized as C(r) ~ , so this corresponds in our notation to H = 1/4, since evidently, 2H — 1/2.) This exponent, 
2H = 1/2, falls well within the experimental range of 0.52 ± 0.12, which is exactly equivalent to the exponent H — 1/4 falling 
within the experimental range 0.26±0.06 which is the one cited in |l^. In the neighborhood of LI, we see that z = 0, indicating 
the relaxation time is independent of length scale and is indicative of a saturated regime. Care must be taken in calculating 
the correlation function since the resulting exponent in (^^ is formally divergent. If we write the correlation function ( |26| ) in 
terms of the auto-correlation function C(t) — 2[S'(0) — S{r)] and use the scaling of the latter as derived in (^) then we see 
that S{r) ~ lim^^o t~^/^ = 0, so that C(r) ~+ 2S'(0) ~ r*^ = const. This reproduces the scaling obtained by Chow and Collins 
by other means. Unfortunately, the error in the null exponent for the saturated regime is unknown. These two power laws 
correspond to the two scaling regions termed as "diffusive" and "saturated", respectively. Note moreover that 11 is repulsive and 
therefore unstable to the slightest perturbation (fluctuation) while LI is stable and attractive. These lines lie within a common 
two-dimensional RG flow domain (in the h — /-plane), and the transition from diffusive to saturated behavior is understood 
from the RG point-of-view as a cross-over phenomenon; see Fig. 6. This plot shows 11 to be at a flnite distance in h above LI. 
The RG flow as determined from the linear model starts from a point on 11 and drops down vertically until it ends up on LI. 
The locus of all possible flows "starting" from points on 11 forms a "curtain" which is depicted in this Figure. We superpose 
the complicated RG flow due to the non-linear interaction for sake of comparison. Actually, while g* — G* — is consistent 
with zero coupling e = and /* — F* holds simultaneously, there is in fact an inverse (and singular) relation between h and 
H, namely, H = 1/h. Since LI has h* = 0, the other line, 11, when plotted in terms of these coordinates {g,h,f), is actually 
infinitely far away from LI This we interpret as an artifact of the model, which remember, is based on an infinitely long string. 
Thus, for an inflnitely long body, the transition to the saturated or bounded phase would never take place (the crossover time 
would be infinite). But this makes perfect physical sense since in fact saturation is a finite-size effect. This connection between 
finite system size and saturation is also drawn for phenomena in surface growth phenomena ^|. Real bodies of course are 
finite in size, and these two fixed lines would be separated by a finite distance in h. 



11 



g 

FIG. 6. Renormalization group flow between 11 and LI in the h — f plane. 



This cross-over time scale is given by a~^, and a fit to the data yields ~ lOsec |L7| . 

Restoring the nonlinear term e > gives rise to additional nontrivial structure in KG parameter space. In total we have 
the two trivial fixed lines {11, LI) plus the four nontrivial fixed points; P2, P3, P4, P5. The scaling behavior of the correlation 
function in both the time and space domains are listed below in Table |[l| . 



Fixed Point 


C(r) = Br2>^/^ 


C{x) = Ax^"^ 


11 


1 




LI 


r» 


X-' 


P2 


^1.5 


x-°-^ 


P3 


^-5.9 




P4 


^-6.0 


x-°-^ 


P5 


^0.04 


x°-' 



TABLE III. 



Scaling behavior of the correlation function about each fixed point. 
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The scaling of C(t) associated with the trivial fixed lines 11 and LI reproduces that obtained previously in where the 
(bare) correlation function was calculated directly in the linear model. The behavior of C(t) in the vicinity of the non-trivial 
fixed points P2, P3, P4 and P5 is due entirely to the presence of the non-linear term ~ y'^ in the equation of motion. The 
scaling behavior in the spatial fluctuations encoded in C{x) was not discussed in but these can be computed just as easily as 
the temporal fluctuations and are listed in the third column of Table 111 for completeness. These together with the behavior of 
C(r) for the non-trivial fixed points constitute predictions of the (non-linear) model within all the phases. We draw particular 
attention to the scaling of the correlation function in the vicinity of the stable spiral P2, which as we demonstrated in the 
previous Section, is an attractor for falling. 



V. RENORMALIZATION GROUP FLOW 



The two-dimensional RG flow for the strictly linear model, Fig 6., is rather featureless and uniform and lends itself to 
easy interpretation. Although we have depicted it in the two-dimensional coupling plane h — f, as already pointed out, 
the phenomenology of posture control picks out a single point on each of the fixed lines 11 and LI, {0,0, F*) and (0,0,/*), 
respectively, and the RG flow is actually a one-dimensional line. By marked contrast, the RG flow for the nonlinear model 
fills out the full three-dimensional coupling space {g,h,f). The relative locations of LI, P2, P3, P4, P5 as shown in Fig. 7, 
where we have suppressed the flow for better visibility. The largest point on the heavy line LI represents the experimentally 
determined point Pec = (0,0,0.05), and one can appreciate the close proximity of the two saddle points P3 and P4, in the 
neighborhood of this point. The stable spiral P2 is further away and the other saddle point P5 is furthest removed from LI. 
The actual numerical coordinates of all these objects are listed in Table |. We have found the RG flow is best represented 
in terms of the (normalized) instantaneous direction field, obtained from directly plotting the differential equations (jl^) as 
functions of {g,h,f). In normalizing, we lose information of the instantaneous "speed" of the flow, but retain the sense of 
flow direction and flow morphology. The instantaneous magnitudes of the flow vectors changes abruptly and dramatically and 
without normalization, makes the graphs extremely difficult to plot. We plot each fixed point with its local system of numerically 
determined eigenvectors. The lengths and directions of each eigenvector system is taken from Table |l| The remaining graphs 
represent a selection of direction flow flelds calculated in the vicinity of the flxed points, for various ranges of the coupling g, 
where the complicated nature of the RG flow can be best appreciated. 




FIG. 7. The fixed line LI, experimental point Pec, and P2, P4, P4, P5. The plotted eigenvectors are taken from Table |. 

The structure of the RG flow is shown in Fig. 8, which gives a "close-up" view of the flow in the vicinity of the experimental 
point and the saddle point P4. Gross feature can be appreciated such as the circulation coming from the region of positive 
/ and sudden change of flow in the region of negative /, which we interpret as signaling the presence of a phase boundary or 
domain. In the closer vicinity of the fixed points the flow is being attracted to the point CC. Since P4 is a saddle point, the 
flow does not end up there, but gets "deflected" when it passes by. 
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FIG. 8. The RG flow in the vicinity of CC and P4. 
In Fig. 9, we show the structure of the RG flow in the vicinity of CC and the stable spiral P2. 



Pec 




FIG. 9. The RG flow in the vicinity of CC and the stable spiral P2. 

In Fig 10, we view the same flow fleld as before, but drawing the points P3 and P4 in their respective positions, and Fig 11. 
gives the view but with extended range in the coupling g, while in Fig 12, the range in g is extended even further and only the 
points Pec and P3 are drawn in place. 




FIG. 10. The RG flow in the neighborhood of CC, P3,P4 and P2. 



14 



g ^P2 

FIG. 11. The RG flow in the neighborhood of CC, P3, P4, and P2. 




FIG. 12. The RG flow in the neighborhood of CC and P3. 




FIG. 13. All points and flow fleld are included. 

VI. SUMMARY AND DISCUSSION 

We have extended Chow and Collin's linear pinned-polymer model of posture control by including a weak quadratic non- 
linearity. There are at least two good reasons for doing so. Firstly, in real anteroposterior movement, the front to back sway is 
not symmetric with respect to the vertical upright position. An obvious way to account for this fact is to introduce a symmetry- 
breaking term in the polymer equation of motion, and the minimal term that can be added is second order or quadratic in the 
displacement field y. Secondly, the effects such as stepping or falling are beyond the scope of the linear model, and these can be 
approximately modelled by means of non-linear terms in the equation of motion . The analysis of non-linear equations is a 
complicated enterprise, but the techniques afforded by the renormalization group permit one to obtain a wealth of information 
regarding the dynamical phases of the system for both large-distance and long-time limits. We have undertaken a detailed 
RG analysis of the fully nonlinear model and have summarized our results in terms of RG fixed points, stability analysis and 
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exponents. We have also numerically computed the fully non-linear RG flow and have represented this in terms of a normalized 
vector flow-field. Knowledge of the RG fixed points is sufficient for determining the exact power-law behavior of the correlation 
function of posture displacement in both the temporal and spatial domains. In the linear limit of the model we recover the 
diffusive and saturated phases of posture control and compute the scaling of the correlation function within each phase. The 
transition from the diffusive to the saturated phase and the associated change in scaling exponent is a crossover phenomena. 
The crossover time is however, finite, only for finite size systems, as we have argued. These results agree with the linear analysis 
of Chow and Collins [Q. The quadratic non-linearity gives rise to four non-trivial fixed points. There are two saddle points 
near the attractive trivial fixed point. While they do alter the RG fiow in the neighborhood of the trivial attractive fixed point, 
they have no bearing on either of the two linear phases (diffusive and saturated) of the model. The linear model fits the posture 
data rather well, and the detailed analysis undertaken here substantiates the claim made in that weak nonlinearities are 
not needed to explain the posture data of healthy individuals. 

The other two non-trivial fixed points consist of an additional saddle point and a stable spiral. The saddle point corresponds 
to a large value of the dimensionless coupling constant and is probably an artifact of perturbation theory. Much more interesting 
is the spiral which is purely attractive and as we have shown, is associated with a falling phase. Its domain of attraction appears 
to be separated from the diffusive and saturated domains. It is important to note that the quadratic nonlinearity is the minimum 
term that can be added to the equation of motion that serves to break the anteroposterior symmetry. This symmetry breaking 
has lead to a falling phase, and thus this one term simultaneously fulfills two distinct requirements. In plf , falling has been 
modelled with a nonlinear potential, but the nonlinearity employed there does not break the y —y symmetry. 

Mention should be made of the short-time inertial effects which are not covered in the present analysis, and which are 
important for an understanding of posture control. These seem to imply the existence of short-time correlations in the noise 
JITI and is no surprise therefore that the purely white uncorrelated noise used here is unable to reproduce this early scaling 
regime. Nevertheless, the analysis carried out here can be straightforwardly extended to handle both white and colored noise, 
and some comments to this effect are provided in Appendix D. 

Although our analysis has centered on the application to posture control, variants of stochastic differential equations of the 
type considered here have applications to a host of other problems where nonlinear waves propagate in a noisy and/or random 
medium |32| , and the general details of the RG analysis carried out here should be useful for addressing these other applications. 
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APPENDIX A: RESPONSE FUNCTION: ONE-LOOP CORRECTION 

The explicit analytic expression for the one-loop correction (hereafter denoted by primes) to the response function ( |l^ ) 
follows immediately from transcribing the diagrammatic representation of the corrected response function (see Fig. 2) into its 
corresponding mathematical elements: 

A'{p,uj) = A{p,iu) + Ae^A{p,Lu) X Ir{p,l^) X A{p,Lj) 

^ A{p,ij){l + Ae^A{p,Lu) X Ir{p,Lj)) 

^ A'-\p,u) = A-\p,uj) ~ 4e^J.(p,a;) + 0(6^), (Al) 

where the loop- integral Ir{p,io) is built up from the bare response function, the bare vertex and the bare noise spectrum 
(for convenience, we have already factored the dependence on the vertex, or bare coupling e, out of the loop integral). From 
inspection of the loop diagram and making use of the Feynman rules, this integral has the structure given by 

Irip,^)^ ^A{p~q,cu-n)C{q,n), 

dq r dn r(g,») 

2-K 2tt (!^2 _^ [^g2 _^ Of _ I3n2]2) - fi] - I3[lj - n]2 + ,y{p - g)2 + a) ' 

valid for an arbitrary Gaussian noise spectral function T{q, Q,). 

The internal momentum and frequency flowing around the loop are denoted by q and fl, respectively. The net momentum 
and frequency flowing into and out of the loop-diagram is p and ui; note that conservation of momentum and frequency is 
maintained independently at each vertex. We take a white noise spectrum r{q, Q.) = 2A, and flrst compute the frequency 
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integral exactly using the residue theorem (the contour may be closed in e ither the upper or lower half-plane). This yields for 
the one- loop corrected inverse response function (from the last line in ( |A1[ )) 

(iuj - 13' u'^ + u'p^ + q') = {iuj- (iuj'^ + up^ + a) + j ^F{q; p, uj) , (A3) 

where the integrand function F, which depends on both internal and external momenta as well as on the external frequency, 
is given by 

F{q-p,uj) = _ I ^^^^^ _ ^^^^^^ _ + _ ^^^^^^ - ^4) } ' ^"^"^^ 

and is expressed in terms of the following poles in the complex frequency plane which arise in the frequency integration: 

^(i- v/4/3(«^9^ + a)-l), 
_L(^_j+^4/3[i.(p-g)2 + a]-l) +uj, 

^(^-i-y/mi4^^WTa\~Tj (A5) 



^2 
^4 



We work with the inverse response function since this is a simple polynomial in p and lo. To renormalize A ^, we must 



expand out the momentum-integral in (A3) in l owest powers in both the external frequency and momentum (a;,p), match 
like-powers on both sides of the expression ( |A3[ ) and then take the hydrodynamic limit tj — > 0, p ^ at the end of the 
calculation. It is important to note that contributions to this asymptotic long-distance and long-time expansion come not 
only from Taylor-expanding the integrand F itself but also from the domain of integration implicit in the integral. We must 
integrate the loop momentum within a fixed "shell" and the net momenta circulating within the loop depends on both external 
and internal momentum variables, and this fact must be taken into account. Thus, the resultant domain of momentum-shell 
integration is given by the intersection of the two intervals A/s < |g| < A and A/s < \p — q\ < A. Up to second order in p, the 
last inequality can be written asA/s+p^ < |g| <A-|-p^, since the 0{p^) terms vanish identically. In taking the intersection 

of this with the first inequality, we have four cases to consider depending on the sign of p (p > 0,p < 0) and the sign of = ±. 
The resultant integration domain, valid for all four cases, can be written as the difference 



dq f dq f dq 



2n I . , 27r , 27r ' 



(A6) 



where the domain E(p, s) = [A/s, A/s -l-p] U [A + p. A] U [A - p. A] U [A/s, A/s - p]. Note of course that E(0, s) = E(p, 1) = (f), 
is just the empty set. 



To proceed with the calculation, in accord with (AS), we need to expand out the function F up to and including quadratic 
powers in both external frequency and momentum {l,uj,p,p'^ ,ujp,uj'^) and consistently combine these with the powers of p 
coming from the integration over E(p, s). In practice, this delicate operation need only be carried out for the renormalization 
of the diffusion constant v. This is because the parameter a does not multiply any positive power of either frequency or 
momentum, so we can take the hydrodynamic limit at the outset in computing its one- loop correction. Next, the parameter /3 
multiplies uj'^, so we must expand the integrand F to this same order to obtain the correction /3' , but we can set the external 
momentum p to zero at the outset: the domain E(p, s) does not depend on external frequency and makes no contribution to the 
renormalization of /3. Terms linear in external frequency (ui) appearing in the loop integral do not yield any new information 
since we can always re-define the time to absorb such corrections when they arise (thus we maintain the unit coefficient 1 in 
front of the term iui in A'). Finally, for the viscosity renormalization, we can set the external frequency to zero at the outset but 
must expand the integrand together with the integration domain up to and including second order in the external momentum 
(p). The constant contribution serves to renormalize a, as we have already remarked. Taking these points into consideration. 



we arrive at the following one- loop expressions for a' ,13' ,v' 
«'^a + 4^£gF(.,0,0) 



'A/s 

e^A dq 1 d'^F{q,Q,Q) 



13^ y^/^ 2n 2! dcj^ 



, 2 2^,e^A/ 2 dq 1 a'f(g,0,0) 



dq dF{q,0,0) 
2n dp 



(A7) 
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Calculating the indicated derivatives of F using the definition of (A4) together with the complex poles (A5) yields the 
following independent equations for the (one-loop) renormalized parameters, namely 

a^2e^A I ' ^ - , 

2n {vq^ + q)2 ' 



A/s 



=P~ e^A 



dq / 
2ti \ 



A/a 
A 



1 



\2(i/i 



27T [AI3{vq^ + q) - 1] Y P 4/3(;.g2 + q) _ 1 - 
3/3^ /3 . 1 



+ 



5/3^ 



{uq^ + a) Ai^vq^ + a)'^ Ai^vq^ + aY 



+ 



(i/g^ + q)2 (z^g^ + a)3 A{vq'^ + a) 



+ 0(S(p,s)). 



(A8) 



These one-loop equations are exact. We have not bothered to explicitly write out the contribution to the viscosity renormal- 
ization coming from the domain S(p, s), since it is easy to show that this will vanish identically when we pass to the differential 
form of the renormalization group equations, i.e, in the limit of a thin-shell. That is, for any function /, we have 



-I 



/(m) du = 0. 



(A9) 



Applying the renormalization group procedure as described in Section 111 (a Kadanoff transformation or coarse-graining 
followed by a re-scaling) to these equations yields the corresponding differential RG equations in (|l7|). 



APPENDIX B: NOISE SPECTRAL FUNCTION: ONE-LOOP CORRECTION 

From the diagrammatic one-loop expansion for the correlation function we obtain the one-loop IPl diagram representing the 
noise spectral function corrected to one-loop, as shown in the diagram in Fig. 4. This translates into the following mathematical 
equation (after factoring out the dependence on the bare vertex or coupling) for a general Gaussian noise spectral function 

V'(p,uj)^V{p,u)+2e^ I^{p,u), (Bl) 

which for the case of white noise considered here, reduces to 

2A' = 2A + 2e^ I„{p,Lo). (B2) 

The loop integral /„ depends in general on external frequency and momentum, whose structure is given by 

Up,^)=r^f ^ C{q,Lu)C{p~q,Q~u;), (B3) 

J J —ac 

but since the noise spectrum is constant, we can take the hydrodynamic limit right away and evaluate the somewhat simpler 
integral 

^"^0' °) - 2n J_ 2. [n^ + [uq^ - + a)2]2 ' ^^^^ 

Once again, the frequency integration can be evaluated easily by the method of residues. Doing so, we obtain the one-loop 
correction to the white noise amplitude 

A' = A + 4e'A' r ^ - r I 3 ^ I 1 (B5) 

J^^^27V A(3{uq2 + a)-l\{uq2 + a) A{uq^ + a)2 A{uq^ + ayJ' ^ ' 

Applying the renormalization group procedure to this equation yields the corresponding RG equation for the noise amplitude 
given in (|l7[). 
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APPENDIX C: VERTEX FUNCTION: ONE-LOOP CORRECTION 



The diagrammatic expansion for the one-loop correction to the vertex function, or couphng constant, is depicted as shown 
in Fig. 3. For general vertex functions, momentum and frequency (i.e, energy) conservation implies that a trilinear vertex can 
depend on at most two independent external momenta and two independent external frequencies. Which two momenta and 
which two frequencies one chooses is immaterial. Translating the vertex diagrams into corresponding mathematical elements 
yields for the one-loop vertex correction the following equation: 

e' = e + 4e^h{ki,LJi-k2,i02). (CI) 

The structure of the one-loop integral is given as follows: 

^ / — {C{q,Q.)A{q~k2,n~u;2)A{q-ki,Q-uJi) 

— OC 

-I- A(-g, -n)C{q -k2,n- LU2)A{q - fci, ^ - (x)i) 

+ A{-q, -n)A{k2 - q,LU2 - n)C{ki - q,LJi - Q.)} . (C2) 

However, since the coupling e is constant, in anticipation of the hydrodynamic limit we can immediately set all external momenta 
and frequencies to zero in computing the one-loop correction to e. Taking this limit, and taking a white noise spectral function, 
we have that the vertex loop integral at zero external momentum and frequency is given by 



7, (0,0; 0,0) = 



^ dq r°° dQ. 2A f \ 



1 



+ 



(in - 130? + vq^ + a){-iQ. - PO? + vq^ + a) 
{-m - (30? + vq^ + afY 

The integral over the internal frequency may be performed exactly, once again by the method of residues. This yields 



e' = e + 8e^y4 [ ^(^—^ , 1 J ^ 3/3 1 

*^ ^ j j^^^2Tx\A{vq^ + a)'-^ 4{uq^ + a)^ {4:l3{iyq^ + a) - 1) \ {uq^ + a) A{uq^ + aY A{uq^ + aY 



(C4) 



Applying the renormalization group procedure to this equation yields the differential RG equation for the coupling listed in 



APPENDIX D: CORRELATED NOISE 



The effect of both temporal and/or spatial correlations in the Gaussian noise spectrum can be also be taken into account in 
this model. Here, we briefly indicate what steps would have to be taken or modifled in the renormalization group program to 
include such correlations. In ^ the spectral function with both uncorrelated (white) and correlated (colored) components is 
written 

V{k,u) = 2A + 2A,.e&-' {^)-\ (Dl) 

where we consider long-range correlations of the power-law type. These are parametrized in terms of two exponents, p and 
9, for spatially and temporally correlated noise, respectively. The naive scaling properties ^ of the stochastic equation are 
extended to include the scaling of the correlated part of the noise, which reads 

A,.o ^ s^^-'-^^+^^^'+''>A,,e. (D2) 

The perturbative expansion for the response, noise and vertex goes through as before, except now, for 6^0, the one-loop 
frequency integrations over Q. must be re-calculated and in general, branch cuts and poles must be dealt with in the complex 
f2-plane. The RGB's will now depend on the two noise exponents p and 6, as will the flxed points and the critical exponents: 
z = z[p, 6), and x = x{Pi There will be an additional RGE for the amplitude of the colored component of the noise yielding 
a total of six equations. By Buckingham's H- Theorem [^, we know that these can be cast in terms of four equations in 
four dimensionless variables. In effect, the correlations in the noise "open up" a new direction in parameter space and yield a 
correspondingly more complicated fixed point and RG flow structure than that of the uncorrelated noise case treated here. 
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